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ABSTRACT 

The long terra evolution of an accreting carbon white dwarf 

( M =* 0»5Mq ) has been studied from the onset of accretion to the 

ignition of helium. The variations in the details of the helium 

shell flash have been examined with respect to variations in mass 

accretion .rate. For intermediate rates ( 10^^-10”^Mgyr~^ ) the 

helium flash is potentially explosive whereas for high rates 
—8 

( 5x10 Mgyr ) the shell flash is relatively weak. The results 
are discussed in the context of the long term evolution of novae. 
Subject headings: stars: accretion - stars; evolution - stars: 
interiors - stars: novae - stars; white dwarfs 



3 


I. INTRODUCTION 

It is generally accepted that the cause of the nova outburst is 
attributable to a thermonuclear runaway in the .hydrogen-rich layers 
of an accreting white dwarf (cf. Starrfield, Sparks, and Truran 
1974a; Prialnik, Shara, and Shaviv 1978). In this model, the hydro- 

i 

gen-rich matter is accumulated until critical temperatures and densi- 

• ••tles *are> reaehed< whereupon «the' accreted ^material 'ds. burned. The 

rapid release of nuclear energy heats the -envelope and eventually leads 

I to its' ejection* The strength of the outburst is primarily determined 
) 

by the abundance of the GNO nuclei and the mass of the white dwarf. 

The long term evolutionary behavior of the white dwarf has recent- 
ly been the subject of several investigations. Two groups of workers 
(Paczynski and Zytkow 1978; Sion, Aciemo, and Tomczyk 1979) have fol- 
lowed a number of hydrogen shell flashes for a wide range of white dwarf 
masses and mass accretion rates. The results of these studies, taken as 
an aggregate, delineate the regime for which hydrogen burning is unsta- 
ble. In this regime, they show that the interflash period is a strongly 
decreasing function of both white dwarf mass and mass accretion rate. 

Work on the underlying helium-rich layer has also received some attention, 
Webbink, Truran, and Gallagher (1978) have proposed that helium runaways 
may be an important ingredient in determining the overall strength of 
the nova outburst. Since the repeated hydrogen shell flashes lead to 
the accumulation of an underlying helium-rich layer, they suggest that 
helium flashes may also be recurrent, but on a longer timescale. They 
point out that the convective mixing which develops during the flash 
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could provide the enhanced CNO abundances in a natural manner (see 

also, Taam and Faulkner 1975). Nomoto and Sugimoto (1977) have also 

focused upon the long term consequences of the helium-rich matter. 

By neglecting the thermal effects of the hydrogen shell flashes on 

the underlying layers , they investigated the development of the helium 

runaway in a helium white dwarf. They found that the character of the 

nuclear flash was extremely sensitive to the efficiency of compressional 

heating and, hence, to the mass accretion rate. For rapid mass accretion 
—8 —1 

( ft i 4x10 «0yr ) , a shell flash resulted which could lead to the mix- 
ing, hypothesized by Webbink et al.- (1978) and for slow mass accretion 

♦ --R —1 

(ft S 2x10 M_yr ) a central flash resulted- leading to a Type I super- 
nova explosion and to the complete disruption of the white dwarf (see al- 
so, Mazurelc 1973).. Although helium white dwarfs in nova binary systems 
cannot be excluded, there is strong observational and theoretical evidence 
which suggests that carbon white dwarfs are more appropriate (see below). 

Accordingly, it is the purpose of this investigation to examine the 
details of the helitma flash in a carbon white dwarf for a range of mass 
accretion rates. In §11, we formulate the problem and outline our basic 
assumptions. The numerical results are presented in §III and in §IV 
we conclude with a discussion of the implications of our results for the 
long term evolution of novae. 

II. FORMULATION AND ASSUMPTIONS 

The accretion of matter is considered in the spherically symmetric 
approximation. Helium-rich material CX“ 0, Y* .98, Z“ .02 ) is ad- 
ded, at a constant rate, to the stellar surface in a quiescent manner 
with an initial specific entropy equal to that of the surface. The con- 
stitutive relations and physical data are similar to those described in 
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Taam' (1979) and will not be presented here. We assume an initial white 
dwarf mass of 0.5Mg, The results of Starrfield, Sparks, and Truran 
(1974b) and Nariai, Nomoto, and Sugimoto (1977) indicate that nova ex- 
plosions are weak or non-existent for white dwarfs of mass less than 
about 0.5Mg. This difficulty .stems from the small surface gravities 
typical of low mass white dwarfs which leads to a low degree of electron 
degeneracy in the hydrogen burning shell and, therefore, to a weak hy- 
drogen shell flash. Since little or no material is ejected from the 
white dwarf surface, a constant mass accretion rate should be applicable. 
Over a long timescale (long compared with the interflash period of the 
hydrogen shell flashes), a significant amount, of helium will, therefore, 
accumulate on the white dwarf. These shell flashes, although important 

g 

for the short term evolution of the nova, are not explicitly included in 
our study. For the mass accretion rates (10”^-10~^_yr**^, see Warner 
1976) which are thought to be relevant to nova binaries, an examination 
of the energy transport from the hydrogen burning shell to the white 
dwarf core reveals that the depth to which energy is conducted is mini- 
mal (Taam 1977). That is, heating of the underlying (helium-rich) layer 
is ineffective as nearly all the nuclear energy is transported to the 
surface. As such, the long term effects of the hydrogen burning shell 
on the thermal structure of the deep envelope layers are, in this case, 
unimportant. On the other hand, the lack of immediate adjustment to 

the newly- accreted material can lead to heating in the interior. For 

-9 “1 

mass accretion rates in excess of about 10 M„yr , coii5>re3sional heat- 
ing will dominate and, hence, the thermal structure will be governed by 
the disparity between the mass accretion and thermal readjustment 
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tdLmescales . The shorter the mass accretion timescale, the more near- 
ly adiabatic will be the compressional heating. In such cases, the 
effects of the hydrogen burning layer may be approximately neglected 
and the problem of the long term evolution of a white dwarf accreting 
hydrogen reduces to one in which helium may be directly accreted. 

As alluded to earlier, a principal conclusion of the study by 
♦ 

Starrfield, Sparks, and Truran (1974a) was the need for enhanced CKO 
abundances in order to reproduce the gross features of fast nova out- 
bursts. This requirement, naturally, led these authors to suggest that 
the outburst occurs on a carbon white dwarf. It is made plausible by 
the fact that the masses of white dwarfs in nova binaries, as. inferred 
from observational data, are large CM'' 1M_, Robinson 1976) coupled 
with the theoretical viewpoint that the cores of the more massive white 
dwarfs should be composed of carbon or material of greater atomic weight. 
The results of single star evolution Cif applicable to binaries) indi- 
cate that stars more massive than 0.45M- ignite helium prior to becom- 
ing white dwarfs (Webbink 1975). Thus, the interior of more massive 
stars should consist of the products of helium burning, i.e. carbon and 
oxygen. It should be noted, however, that more massive helium -white 
dwarfs (M > 0.45Mg) could be built up by mass accretion in binaries-. 
Normally, the production of a white dwarf requires a previous giant 
stage and, therefore, a relatively wide initial separation between 
the components of the original binary system. Since nova binaries 
are found at short periods ( "several hours), the manner in which such 
systems are formed remains obscure. Evolutionary processes by which this 
transformation may have taken place are discussed by Ritter (1976) and 
Faczynski (1976). Since the details of ,the mechanisms have yet to 
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be quantified, we adopt the point of view that we may investigate the 
evolution of a white dwarf in a short period system without regard to 
the issue of its formation. 

III. RESULTS 

The evolution of a carbon white dwarf accreting helium-rich matter 

i 

was followed utilizing the stellar evolution code developed and described by 
by Eggleton (1971, 1972). This hydrostatic code utilizes a non-Lagranglan 
mesh, which includes the entire envelope, ■ and is wall suited for investiga- 
tions of mass accretion (or mass loss) . The initial model is characterized 
by an effective temperature T = 1.23xlo\, a luminosity L = 4.31xl0“\„, 
and a central temperature 'T^ — 1.8x10 K. Three evolutionary sequences 
were computed, each distinguished by the rate of mass accretion (M ='.l, 1, 
5x10 ^M^yr ^) . . 

For these rates-, compressional heating effects determine the course 
of the evolution. As mass is accreted, the inner layers contract and 
the gas is compressed, thereby releasing gravitational energy in the pro- 
cess. Since the mass accretion timescales are shorter than the thermal 
readjustment timescales in the interior and the deep envelope layers, 
this energy cannot be radiated quickly enough and the interior tempera- 
tures rise. In fact, compressional effects significantly affect the 
central regions of the white dwarf as well. This is in contrast -to the 
static interior conditions found in the short term evolution because 
the amount of mass accumulated (AM > O.lMg), here, is not negligible, 
conqjared to the total mass of the white dwarf. In the following, we 
present the results for the cases of rapid, intermediate, and slow mass 


accretion.' 
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(1) Rapid Mass Accretion 

—8 ••1 

A high rate of mass accretion ( fl = 5x10 M„yr ) is considered 
in order to maximize the effects of compressional heating. The numer- 

I 

ical results of the evolutionary calculations are summarized in Table 
1 and illustrated in the temperature-density diagram in Fig. 1. It is 
clear that the thermal structure of the entire white dwarf is affected. 

In particular, an off-centered temperature maximum forms in the vicinity 
of the He-C composition interface (corresponding to the region of par- 
tial degeneracy) where the gravitational energy release associated with 
•the compression is effectively localized. Note that the mass point cor- 
responding to the temperature maximum does not coincide with the base 
of the accreted layers (see Fig. 1). For, the high, degree of electron 
degeneracy C i 100, where is the electron degeneracy parameter) char- 
acteristic of the interior, the high efficiency of electron conduction 
inhibits the growth of large temperature differences. The more effec- 
tive heating in the partially degenerate regions is made^^dent upon 
e x a mi nation of Fig. 1 where it is seen that the rise in shell tempera- 
ture, T^, increases from S.ixlO^K to I.lxlO^K in 7.9xl0^yrs. as com- 
pared to the rise in central temperature, T^, from 1.9xl0^K to 2.Xx 

lO^K (see Table 1). After 0.112Mg has been accumulated ( t ~ 2.24xl0^yrs) , 

5 —3 

the physical conditions are such that (p = 1.7x10 gm cm and T = 8.8 x 

m m 

10 K) helium is Ignited. Because of the strong temperature dependence of 
the helium energy generation rates and the moderate degeneracy of the 
electron gas in the burning shell (t|) " 7) , a degenerate shell flash en- 
sues* Within 1.6x10^ yrs., a convective region forms immediately above 
the nuclear burning shell. At this stage, the shell is characterized 
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8 S ^ 

by a temperature of 1.1x10 K, a density of 2.05x10^ gm and a nu~ 

3 

clear lu min osity of 2.9x10 Lq. It takes an additional 40 yrs. to reach 
8 8 

1.8x10 K and 6.9x10 Lq by -which point the degeneracy of the electron 

gas is decreasing. At this time, convection is inefficient as super- 

adiabatic temperature gradients were required to transport the flux. 

The mass in the growing convective shell increased to 0.08M„, and the 
♦ - 6 - 

3 —3 

densities at its outer boundary decreased to 4x10 gm cm . The subse- 
quent evolution of the helium burning shell is, now, governed by the 
^ nondegenerate shell flash instability. In the next 5.6x10^ seconds, 

1 the density, temperature, and nuclear luminosity in the shell reached 
10 gm cm , 2.5x10 K, and 7.3x10 Lg, respectively. As a result of the 
expansion of the outer core (i.e. He-C composition boundary), the central 
regions expanded and cooled adiabatically' (see Table 1). At all phases 
of the thermal flash, the star remained in hydrostatic equilibrium. The 
calculations were discontinued when it became apparent that the progress 
of the shell flash was qualitatively similar. to that fomid“in case A of 
•Nomoto and Sugimoto (1977). Although Nomoto and Sugimoto (1977) con- 
sidered a 0 . 4Mg. helium white dwarf , a .comparison of Table 1 with their 
Tables 2 and 3 shows reasonably good agreement. 

(ii) Intermediate Mass Accretion 

— 1 

For this evolutionary sequence, a mass accretion rate of 10 M..yr~ 

0 

is considered. Because of the longer mass accretion timescale, the com— 
presslonal heating is less adiabatic, and the energy is not as effective- 
ly localized as found in case (i) . The results of our computations to 
the ignition of hellm are presented in Table 2. Altlwugh .the compres- 
sion rate considered, here, is lower, the central temperatures exceed 
those found to develop in case (i) . This is attributable to the greater 
amount of mass which is accumulated prior to the Ignition of helium 
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(0.617M as compared with 0.12M for case (i) ). It is clear from 
o 0 

the temperature-density diagram (Fig. 2) and from Table 2 that after 
0.55 Mq has been accreted, the heat transport in the interior results 

in an isothermal distribution. Helium ignition occurs at a density of 

7-3 7 

p =» 1.45x10 gm cm and a temperature of T = 7.35x10 K. Because of 
m m 

the high degree of electron degeneracy (i|; ' 100) , this situation be- 
comes potentially explosive. Since the pressure is insensitive to var- 
iations in the temperature, the nuclear burning will occur under hydro-- 

static conditions until the heating timescale, t , is reduced to val- 

nuc 

ues comparable to the hydrodynamic response timescale, t -- (see below) . 

rf 

That is, under these circumstances, stellar expansion and adiabatic cool- 
ing cannot control the nuclear reactions (cf. Arnett 1969). As a result, 

the shell temperatures rise rapidly. For instance, the maximum shell 

8 8 

tenq)eratures increased from 3.3x10 K to 6.24x10 K in 0.002 seconds as 

8 8 

compared to the rise in temperature from 1.21x10 K to 1.57x10 K in 6.55 

3 7-3 

xlO seconds. Even at the high densities (p ' 10 gm cm ) characteristic 

of the nuclear burning region, the convection could no longer be con- 
sidered efficient. Superadiabatic temperature gradients (7 - V ) of 

ad 

about 4 were required to transport the energy. These steep temperature 
gradients necessarily confined the rapid burning to a narrow shell (the 
energy generation rates decreased, by more than 6 orders of magnitude for 
a 30% decrease in pressure) . It can be seen from Table 2' that (i) the 
characteristic timescale for the temperature rise in the burning shell 
due to nuclear reactions, denoted as becomes increasingly short 

with greater temperatures and (ii) the expansion of the shell is small 
(ALogp^ •* 0.156). A con 5 )arison of with the free fall time- 

scale which we take to be approximately equal to the adiabatic expansion 
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timescale, at log T = 8,795 yields t / =* 0.01. Under < 3 uch con- 

ditions, significant burning is possible prior to the hydrodynamic res- 
ponse which must ultimately govern the further evolution of the star. 

1 

At this point, the calculation was discontinued as hydrodynamic 
effects, which were not included, were starting to become impor- 
tant. The thermonuclear flash develops on such a short time scale 
. .that .little energy is. conducted to. the central regions. It is 
expected that during the hydrodynamic phased expansion and adiabatic 
cooling of the central regions will result. For the physical condi- 
tions in the helium burning shell found here, it is possible that a 
detonation wave will form offcenter and propagate through the helium- 
rich material. This is in contrast to the investigations of Mazurek 
(1973) and Nomoto and Sugimoto (1977) where a detonation wave formed 
at the center of a helium white dwarf , 

(ill) Slow Mass Accretion 

In order .to exaggerate the differences between the results of se- 

—9 —1 

quences (i) and (ii) , a mass accretion rate of 10 M^yr was chosen. 

The evolutionary changes in the structure. are exhibited in the tempera- 
ture-density diagram (Fig. 3) for the phases prior to helium ignition. 
During the initial stages of accretion, it is seen that the temperature 
distribution becomes non-mono tonic, but after this transient phase, the 
white dwarf evolves without a temperature inversion. Because the accre- 
tion timescale is 10 times longer than the intermediate case, the com- 
presslonal heating is not as effective and the onset of the degenerate 
flash is delayed to much greater densities (p " 3x10 gm cm , 


see Table 3) . 
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Because of the strong electron screening factor, the helium burning 

8 7 

rates were enhanced by a factor of 10 C at T =* 5.3x10^ K) . Since 
the densities in the shell are. more than an order of magnitude greater 
than found in case (ii) , the time scale for the development of the 


flash-, from a given temperature, would be greatly reduced. Helium 


ignitiotf was found to occur after 0.88 M had been accumulated. The 

■ ■ ® , ■ 

calculations were followed to the onset of shell helium burning 


and not through the degenerate shell flash' stage as the thermal 
jevolution is expected to be similar in kind to that calculated 


for case (ii) . 


IV. CONCLUSION 

Thermonuclear runaways initiated by the accretion of helium-rich 
matter onto a carbon white dwarf of 0.5 Mg have been investigated to ap- 
proximate the long term evolution of novae. For the range of mass ac- 
cretion rates studied ( M = 0.1-5x10 ® M yr compressional heating ef- 
» 0 

fects were found to be significant in raising the interior temperatures 

to values appropriate for the ignition of helitim. The effectiveness of 

this heating dictated the subsequent behavior of the helium shell flash. 

Generally#- the lower the rate of compression, the greater the strength 

of the shell flash as helium ignition is delayed until a greater density 

(electron degeneracy) is reached. For M = 5x10 ^ M„yr- ^ the flash pro- 

. -8 -1 

ceeds under hydrostatic conditions, whereas for M** O.l-lxlO ^ 0 ^^ 

the helium may detonate. In the former case, the extensive convection 

zone formed exterior to the main energy generating region extended to 

3 -3 

densities ( p- 4x10 gm cm ) which may lead to the mixing of carbon 
into the overlying hydrogen rich layers (see also, Webbink et al. 1978). 
This mechanism could provide the CNO enhancement required for fast 
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nova outbursts (cf. Starrfield, Sparks, and Truran 1974a). The hydro- 
gen shell flash immediately following the helium runaway is expected 
to be most rapid and the strength of each succeeding hydrogen flash to 
be less intense as the abundances of the CNO nuclei are reduced. A 
similar progression of hydrogen flashes would be expected after each 
-'Succeeding -helixim runaway • ( t ~ ’-‘2r4xl0^ -yrs) . '“"For^'H = 10 ^ 
the complete combustion of all the helium matter is likely. The en- 
ergy yield would be of the order of 10^^ ergs which is greater than 
the binding energy of the helium matter by about a factor of 4. If 
the detonation wave forms and propagates to the surface, the result 

could be a supernova explosion. Since the time scale for the evolu- 

4 

tion is short (tS 10 sec) little heat is conducted into the central 
regions, and it is quite .possible that the carbon white dwarf core will 
not participate in the .evolution, and therefore, remain intact. For 
the lowest mass accretion rate studied ( M » 10 )» likli- 

hood of the formation of a detonation wave is reduced as the internal 
energy of the degenerate electrons at the burning shell is comparable 
to the nuclear energy release. Only a fraction of the envelope may 
be burnt » perhaps resulting in a less energetic event (a nova). Since 
compressional heating is much less adiabatic for lower mass accretion 
rates ( M ^ 10 ) , helium burning in the white dwarf interior 

may be avoided in such cases. On the other hand , energy transport from 
the hydrogen burning shell is non-negligible and it is possible that 
helium burning initiated by the hydrogen shell flash will develop in 
the outer surface layers (Taam and Faulkner 1975, Team 1979). Investiga' 
tions along these lines are clearly needed. 
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Finally we point out that mass ejection during the hydrogen shell 
flash phase occurs for high mass white dwarfs even with normal- CNO 
abundances (Sparks, Starrfield. and Truran 1978, Prialnik, Shara, and 
Shaviv 1978). Here, radiation pressure plays the major role in ejecting 
the matter. Our results, which have been derived assuming a constant 
mass acci^etion rate,. should still be applicable as, long as a time aver- 
aged rate is- used. 


This research has been supported by NASA through grant NGL 14-007-41 
and by Northwestern University. The numerical computations were performed 
at the computing center of Northwestern University. 



TABLE 1 


EVOLUTION AT H - SKlO'^M^yr 


t (years) 

4.72(5) 

1.27(6) 

2.24(6) 

2.40(6) 

2.40(6) 

2.40(6) 

2.40(6) 

log 

4.608 

4.981 

5.239 

5.326 

5.288 

5.221 

5.003 


7.713 

7.853 

7.942 

7.985 

8.140 

8.230 

8.404 

AM(Mq) 

2. 36 (-2) 

6. 33 (-2) 

1. 12(^1) 

1. 20(^1) 

1.20(-1, 

1.20(-1) 

1.20(-1) 


4.03 

5.22 

6,28 

.48 

4.15 

'2.89 

0.84 

log 

6.281 

6.364 

6.465 

6.480 

6.476 

6.462 

6.435, 

log T^ 

7.284. 

.7.332. 

7.393 

7.402 

7.400 

. 7.391 

7.374 

♦ 

T (s) 

nuc 

- 

- 

4.00(14) 

3.76(12) 

4.57(7) 

4.36(5) 

1.56(3) 

e (erg. 


_ 

1.83(1) 

2.14(3) 

2.52(8) 

3,25(10) 

1.36(13) 


T^ represents the maximum temperature in the helium-tlch region, 

p represents the density at the mass point corresponding to T , and 
m . . m 

represents the maximum energy generation rate. 



TABLE 2 


EVOLUTION AT M “ 10“®MQyr”^ 


t(years) 

9.79(6) 

2.09(7) 

4.26(7) 

5.50(7) 

6.17(7) 

6.17(7) 

6.17(7) 

.6.17(7) 

6.17(7) 

log 

5.170 

5.600 

5.933 

6.207 

7.161 

7.102 

7.077 

7.059 

7.005 


7.703 

7.804 

7.807 

7.815 

7.866 

8.024 

^321 

8.518 

8.795 

AM(Mq) 

9. 79 (-2) 

2.09(-l) 

4.26(-l) 

5.50(-l) 

6.17(-1) 

6.17(-1) 

6.17(-i) 

6.17(-1) 

6.17(-1) 


7.88 

14.88 

23.69 

33.75 

96.15 , 

62.71 

30.70 

19.02 

.9.31 

log 

6.443 

6.685 

7.186 

7.502 

7.689 

7.689 

7.689 

7.689 

7.689 

log T^ 

7.439 

7.723 

7.783 

7.825 

7.856 

7.856- 

7.856 

7.856 

7.856 

T (s) 

nuc 

.. 

- 

- 

- 

4.53(13) 

1.53(7) 

1.82(1) 

1.25 (^2) 

7.67(-4) 

e„(erg g“^a”^) 

•• 

••• 



1.35(2) 

5.76(8) 

9.57(14) 

2.20(18) 

6.76(19) 



TABLE 3 





EVOLUTION AT 

M = 10"^Mgyr‘ 

-1 



t (years) 

9.44(7) 

2.24(8) 

4.27(8) 

6.12(8) 

7.23(8; 

8.17(8) 

8.77(8) 

log P„ 

5.446 

5.983 

6.505 

7.186 

7.526 

7.997 

8.478 

■'m 

7.414 

7.489 

•7.512 

7.546 

7.583 

7.635 

7.721 

AM(Mq) 

9. 44 (-2) 

2.24(-l) 

4.27(-l) 

6.12(-1) 

7.23(-l) 

8.17(-1) 

8.77(-l) 

'i' 

2.94(1) 

5.29(1) 

9,96(1) 

2.07(2) 

2,74(2) 

3.90(2) 

5.02(2) 

log 

6.447 

6.739 

7.206 

7.734. 

8.115 

8.501 

8.899 

log 

7.393 

7.505 

7.525 

7.558 

7.596 

7.648 

7.723 

T (a) 

nuc ' 

- 


- 

- 

- 

- 

1.35(13) 


EgjCerg g'^s S 


3.23(2) 
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Figure Captions 

Figure 1. Structural changes in the temperature density plane from 

• —8 

the onset of mass accretion to the ignition of helium for M = 5x10 M^yr 
The time corresponding to each curve is given for clarity. The uniabel- 
led curve corresponds to the distribution in the initial white dwarf 
model, " Vertical lines Indicate ‘ the" location of the He-G composition 
interface. 

—8 —1 

Figure 2. Same as for Figure 1 except M =» 10 Mgyr . 

-9 -1 

Figure 3, Same as for Figure 1 except M = 10 Mgyx . 
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